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We demonstrate computerized tomography (CT) reconstructions from absorption, phase and dark- 
field signals obtained from scans acquired when the x-ray probe light is modulated with speckle. 
Two different interlaced schemes are proposed to reduce the number of sample exposures. First, the 
already demonstrated x-ray speckle-vector tracking (XSVT) concept for projection imaging allows 
the three signal CT reconstructions from multiple images per projection. Second, a modified XSVT 
approach is shown to provide absorption and phase reconstructions, this time from a single image 
per angular projection. Reconstructions from data obtained at a synchrotron facility emphasize the 
potential of the approaches for the imaging of complex samples. 
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I. INTRODUCTION 

X-ray computed tomography pQ has become an invalu¬ 
able tool for nondestructive testing and 3D rendering of 
samples in fields as broad as material sciences, cultural 
heritage, paleontology and even medical science. In the 
latter domain, the dose necessary for imaging in vivo 
samples often becomes a critical issue and so must be 
carefully controlled and maintained below an acceptable 
hazard threshold m- Henceforth, phase and dark-field 
contrast-imaging methods are regarded as sensible ap¬ 
proaches to eventually enhance the rendering quality and 
quantity of information on soft tissues with a dose level 
equivalent to or lower than that encountered by using 
conventional absorption contrast radiography. 

Nowadays, a few methods and instruments are avail¬ 
able to recover the phase shift induced by a sample ex¬ 
posed to an x-ray light wave. Most techniques use the 
partial transverse coherence property of an x-ray beam 
to generate interference between optical waves and en¬ 
able the calculation of the phase shift. This applies, for 
instance, to propagation-based methods relying on either 
the contrast transfer function (holotomography) [4] or on 
the transport-of-intensity equation [5]. In parallel, the 
demonstration of grating-based phase-contrast imaging 
using a laboratory x-ray source promotes the attractive¬ 
ness of this method [6]. Therein, a clever setup enables 
an extended source to be split into a series of smaller- 
sized sources to permit the generation of constructive in¬ 
terference compatible with shearing interferometry 00 . 
Another method sensitive to the differential phase of an 
x-ray beam is analyzer-based imaging that employs crys¬ 
tals with a narrow Darwin bandwidth to render different 
contrast for rays having different propagation orienta¬ 
tions [9]. 

Near-field speckle-based methods Q2HH form a more 
recent class of x-ray imaging techniques and are proven 
to be applicable both at synchrotrons and at laboratory 
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sources na EE]. Their attractiveness lies in their low 
requirements on coherence and in the simplicity of the 
wave-front modulator, reduced to a simple diffuser such 
as a piece of sandpaper or a biological filter. 

Three main speckle methods relying on pattern cor¬ 
relation in the near field are available today. We pur¬ 
posely leave aside the near-field ptychography technique 
m that is based on a different principle (iterative phase 
recovery) but which often employs a scattering object to 
improve its efficiency. The first introduced method is x- 
ray speckle tracking (XST) [llj 12. [15] ; therein, small 
subsets of speckle are tracked between two images, one 
being taken with the sample present in the beam, and the 
other being used as a reference. The second method re¬ 
quires a slightly more complex arrangement, as the scat¬ 
tering membrane needs to be scanned transversally to 
the beam with micrometer- or submicrometer-scale steps. 
This x-ray speckle-scanning (XSS) method [13] provides 
very high sensitivity at the cost of a number of sample 
exposures that can quickly rise up to a few dozen. In the 
last and most recent method m, the scattering object 
used to generate the static speckle is moved to various 
transversal positions with respect to the beam while vec¬ 
tors built from speckle images are tracked across refer¬ 
ence vectors. We decide to label this method the x-ray 
speckle vector tracking (XSVT) technique on account of 
the aspects it shares with the XST technique and the 
vectors involved in it. 

In the first section of this paper, we briefly review 
and compare these three speckle-processing approaches. 
Efforts are made to highlight their differences and to 
demonstrate theoretically, and later experimentally, that 
XSVT is particularly well suited to performing 3D mul¬ 
timodal computed tomography. Indeed, without any a 
priori assumption on the samples, only a few images per 
projection are required for multimodal tomography when 
employing the appropriate scheme. In contrast to previ¬ 
ous work nn no], the XSVT approach provides 3D ren¬ 
dering of the dark-field signal, in addition to the absorp¬ 
tion and phase signals, whilst the spatial resolution of 
the imagery is limited only by the detector used. 

In the last part, we introduce a mixed approach based 
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on the combination of the XST and XSVT methods. This 
original scheme permits the extraction of the phase sig¬ 
nal with a resolution approaching that of the XSVT tech¬ 
nique with a lesser number of exposures than for XSVT 
alone. Experimentally, the phase-contrast tomography 
reconstruction obtained from this mixed approach shows 
that the number of sample exposures necessary is com¬ 
parable to the case of absorption tomography. 


II. SPECKLE-PROCESSING SCHEMES 

We give here a short comparative review of the speckle¬ 
processing approaches. 

The three speckle schemes available today differ in the 
kind of data collected, the numerical treatment applied 
and, for two of them, their intrinsic sensed signals. Nev¬ 
ertheless, the concept shared by all three schemes is the 
use of cross-correlation operations to track the displace¬ 
ment of a distorted pattern with respect to a reference 
pattern [21] caused by x-ray refraction through the sam¬ 
ple. In all schemes, the speckle pattern used as struc¬ 
tured illumination is generated by interference upon the 
propagation of coherent photons through a scattering ob¬ 
ject (scatterer) with small random features, often made 
of sandpaper or of a biological filtering membrane with 
a known statistical feature size. 

To be more precise, the insertion of a sample into an 
x-ray beam generates a dephasing of the beam waves by 
the amount </> = —k f S(x, y, z)dz where k is the wave 
number and S is the optical refractive index of the sam¬ 
ple material. Upon propagation along an e z axis, this 
dephasing is responsible for the distortion of an input 
reference modulation pattern, i.e., here, the speckle pat¬ 
tern. The local speckle shift is equal to v = in 

the basis (e x ,e y ) transverse to the photon-propagation 
direction. Therein, V is the del operator and d the prop¬ 
agation distance from the sample to the detector. Note 
that this equation is valid in the x-ray regime due to the 
small diffraction and refraction angles involved with such 
light. 

This linear relationship between the local displacement 
v of the modulation features by comparison to a refer¬ 
ence position on one side and the phase gradient induced 
by the sample on the other provides a simple and com¬ 
mon way of accessing beam phase information. This core 
principle is, for instance, applied within the grating inter¬ 
ferometer [7 and coded aperture [22] approaches. Here, 
the use of near-field speckle as a wavefront modulator is 
particularly beneficial in the hard x-ray regime because 
such a pattern is shown to be invariable upon propaga¬ 
tion over distances that can range from millimeters to 
meters m- Furthermore, this kind of speckle is easily 
accessible even with beams presenting low transverse co¬ 
herence properties. 

The distinctive aspects of each speckle scheme lie in the 
approach employed for the recovery of v(x,y). Figure [l] 
illustrates the data that are collected and used within 


each method. 

Within XST (Fig. [I] (a)) [HIHI], the scattering object 
is fixed and only two images are necessary: one with and 
one without the sample present in the beam path. The 
processing consists of taking a small subset of pixels con¬ 
taining some speckle pattern and tracking its shift from 
the sample image position across the reference speckle 
image using a cross-correlation algorithm. The opera¬ 
tion is repeated for all subsets centered on each pixel. 
Thus, the vector map v(x,y) is calculated with subpixel 
accuracy and the two orthogonal differential wave-front 
gradients are recovered. Generally, the size of the tracked 
subset can range from 5x5 pixels to more than 21 x 21 
pixels. Here the resolution is limited by the grain size of 
the scatterer, as the superimposition of the speckle grains 
onto sample features of equivalent size renders the track¬ 
ing of the modified speckle pattern unstable. Despite 
the existence of algorithms that take into consideration 
the effect of subset distortion to enhance the tracking ro¬ 
bustness [18], the method gets, locally, less robust and/or 
accurate in the presence of a strong phase shift, e.g. at 
the sample boundaries. 

With XSS [13], whose principle is shown in Fig. [TJ (b), 
the scattering object is mounted on a highly accurate 
bidirectional translation stage in order to raster scan the 
scatterer transversally to the x-ray propagation direction. 
Two sets of images are then collected by performing the 
same scan, whose step is on the order of 100 nm to a 
few micrometers, with the sample present in the beam 
and without it (reference). For each pixel, as marked by 
an eye symbol in Fig. [l] (b), two data arrays are col¬ 
lected and the analysis is done pixel by pixel. A cross¬ 
correlation-based treatment computes the shift v{x,y) 
between the recorded patterns with a substep accuracy. 
Finally, the differential phase gradient is retrieved after 
linear normalization by a geometrical factor [14]. This 
scheme can offer a very high spatial resolution when com¬ 
bined with a strongly magnifying setup. Moreover, since 
the definition of the wave-front gradient is here propor¬ 
tional to the step size employed, it is possible to reach 
single nanoradian accuracy. On the other hand, such per¬ 
formance is accessible only at the cost of many sample 
exposures due to the necessity of performing 2D mesh 
scans, or at least, two wide-range ID scans [T3]h 

XSVT is the most recently introduced speckle ap¬ 
proach M- It requires the mounting of the scattering ob¬ 
ject on a moving stage capable of high repeatability with 
respect to the effective pixel size of the detector upon geo¬ 
metrical demagnification brought by the setup. However, 
unlike XSS, no accurate calibration on the motor step size 
is required. Figure [T] (c) shows the stacking up of images 
when the sample is present and when it is removed from 
the beam. The stage repeatability is needed for the scat¬ 
terer to be in exactly the same position when recording 
images with and without a sample. Under this condition, 
the data-acquisition procedure consists of recording each 
r] th pair of sample and reference images of the two stacks 
with the scattering object located at every given posi- 
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FIG. 1. Conceptual sketches showing the data recorded and processed within each speckle-processing scheme, (a) X-ray Speckle 
Tracking, (b) X-ray Speckle Scanning, (c) X-ray Speckle-Vector Tracking. The reference data are drawn in red (top line) and 
the ’sample’ data in blue (bottom line). 



XST 

XSS 

XSVT 

Specific setup 
requirement 

None 

Accurate 

motor 

Repeatable 

motor 

Sample exposure 

1 

o(v 2 ) 

o(v) 

Sensitivity 

~ 0.1 /xrad 

~ 0.001 /xrad 

~ 0.01 /xrad 

Spatial 

resolution 

~ 10 fim 
( grain size) 

~ 0.1 /im 
( pixel size) 

~ 1 ynri 
(> pixel size) 


TABLE I. Comparative table of speckle based techniques. 


tion of the scatterer. It is also necessary that a fixed pair 
of images possesses locally a speckle pattern that is un¬ 
correlated from those of the other images. Data analysis 
starts by building two speckle vectors for each pixel using 
the two image stacks. Then, each vector from the sample 
stack is tracked across those of the reference stack. The 
location of the peak of maximum correlation between the 
samples tracked vectors and the reference ones provides a 
mean of recovering the local displacement vectors v(x,y) 
and hence the differential phase gradient induced by the 
sample. This protocol means that a resolution close to 
that of the detector can be reached and that the sensitiv¬ 
ity scales with the detector pixel size and the propagation 
distance. 

Table [I] presents a brief survey of the attributes of each 
of the three techniques. Note that the figures provided 


are discussed in more detail elsewhere hh nans na 123 - 
El]. We recognize that, in all three methods, the dark- 
field signal can be obtained by considering the change in 
speckle visibility between the references and sample im¬ 
ages. However, although this capability is demonstrated 
with XST [15], only the XSS and XSVT approaches 
presently prove to be quantitative and accurate enough 
for 3D tomographic reconstruction. The moderate num¬ 
ber of images per projection necessary for XSVT to access 
the 2D differential phase gradient in addition to the dark- 
field signal constitutes the main asset of the technique. 
The data redundancy in the XSVT technique permits the 
sensitivity to be improved down to 10 nr ad. From these 
considerations and from Table [T] we are able to conclude 
that the XSVT method represents the best compromise 
on the number of necessary images per projection with re¬ 
gard to the sensitivity and spatial resolution. Therefore, 
although tomography is demonstrated with XST mm, 
XSVT represents our approach of choice to achieve high- 
resolution phase-contrast tomography at a cost-effective 
number of sample exposures. Nevertheless, one can fore¬ 
see that XSS will be able, in the future, to perform 3D 
tomographic imaging of samples presenting very small 
electronic variations when the dose is not an issue. 
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III. EXPERIMENTAL SETUP 



or a (I)/ < I >= 0.13. The recorded images have a size 
of 2048 pixels horizontally by 1850 pixels vertically. 

The selected samples presented here are juniper berries 
and red currant berries. The juniper sample is interesting 
for the presence of features at very different spatial fre¬ 
quencies and for the significant amount of x-ray scatter¬ 
ing generated by structures at the submicron scale. By its 
very nature the juniper berry is a difficult sample to im¬ 
age for many techniques. In comparison the red currant 
presents fewer parts that generate strong x-ray scattering 
but has large homogeneous features that can sometimes 
become problematic as they render less-accurate quanti¬ 
tative results with propagation-based methods. 


IV. XSVT-BASED TOMOGRAPHY WITH AN 
INTERLACED SCHEME 



FIG. 2. Bottom: setup of the experimental arrangement. 
Top: acquisition scheme used for the data collection. For 
each angular projection, a set of four images images (either 
A, B, C, D or E) are collected with the membrane located at 
distinctive positions. Each of these sets of 4 different speckle 
illuminations are reused every five projection angles, meaning 
that in total 20 independent speckle illumination are used. 

A sketch of the experimental setup is shown in Fig [2] 
The experiment is conducted at the ESRF beamline 
BM05 [26], where a precision tomography station is per¬ 
manently installed. At BM05, the photons are produced 
by a 0.85 T bending magnet on the circulating 6.02 GeV 
electrons of the storage ring. The photon energy is nar¬ 
rowed down to a band of A E/E ~ 10 _4 centered around 
17 keV by means of a double-crystal Si(lll) monochro¬ 
mator. The sample is placed at 55 m from the source and 
the detector, a CCD-based FReLoN (Fast Read-out Low 
Noise) camera coupled to a scintillator, d=lm further 
downstream from the sample. The effective pixel size of 
the imaging system is h = 5.8 pm. A piece of abrasive 
paper with a grit designation of P800 is mounted on a lin¬ 
ear translation stage allowing displacements in the beam 
transverse directions and located d p = 400 mm upstream 
of the sample. At a distance of d = 1 m from the sample, 
denoting by< > and cr, respectively, the mean value 
and the standard deviation of the vector scalar compo¬ 
nents, the visibility of the speckle presenting a nearly 
Gaussian distribution is of {Imax dmin )/2 < I >— 0.4, 


The method of structured illumination employed in 
this section has been quickly presented in Sec. [TT] and 
is detailed for projection imaging Ref. m- Here in addi¬ 
tion, an interlaced scheme is used to optimize the number 
of images collected for each projection. The proposed in¬ 
terlaced XSVT method allows us to greatly reduce the 
overall number of sample exposures whilst conserving the 
accuracy of the phase and dark-field signals by factoring 
in the images recorded at consecutive projection angles. 

The data-acquisition scheme is sketched at the top of 
Fig-U a total of 20 distinctive transverse positions of the 
membrane with respect to the beam direction are defined 
and split in five sets p = [A, F>, C, D, E\. For each of the 
N = 1800 angular projections over a 180-degree tomog¬ 
raphy scan, four images are recorded with the membrane 
located at the positions k p = 1, 2, 3,4 defined for each an¬ 
gular set. For each projection angle Q n with n £ | [1,7V]|, 
these positions are defined as the modulo 5 of the pro¬ 
jection number thus generating a cyclic repetition of the 
set p employed. 

Experimentally, the full data set results from the 
collection of 20 different scans of 360 projections, the 
starting scan angle being shifted by 0.1 degree every 
four scans. Then, for each projection n, the interlaced 
scheme consists of using the image intensity i(fl n ,x,y,k p ) 
recorded at the projection {n — 2,n — l,n, n + l,n + 2} as 
a multivariate random variable to build speckle vectors 
sof size 20: 

s v (f4 n , x, £/) = ^^i(U n _ 2 , X, y, kp) ^ 

\ ^(U n , x, y, k p ) \ ,..., (1) 

^ ) k p 

i(U n _|_2} x, £/, kp) i 

) kp 

with the sets denoted { } built with k p = 1, ..,4. 

Such an interlaced scheme ensures the avoidance of 
any redundancy in the speckle vectors that would other¬ 
wise decrease the amount of independent statistics they 



















5 



0.5 1 1.5 2 2.5 0.5 1 1.5 -10 0 10 -10 0 10 


Raw image (arb. units) Transmission (arb. units) Horiz - Wavefront grad, (grad) Vert. Wavefront grad, (grad) 



-10 0 10 -10 0 10 -600 -400 -200 0 0.6 0.8 1 

Horiz. Wavefront grad, (grad) Vert. Wavefront grad, (grad) Phase (rad) Scattering (arb. units) 


FIG. 3. Images of the juniper-berry sample, (a) A raw image of the juniper-berry sample at the first angular projection with an 
inset showing the speckle modulation of the wave-front, (b) The corresponding calculated absorption image, (c) Horizontal and 
(d) vertical wave-front gradients calculated with the XST technique [T T\. (e) Horizontal and (f) vertical wave-front gradients 
calculated with the interlaced XSVT method, (g) Phase image from the 2D integration of (e-f). (h) Dark-field image obtained 
with the XSVT interlaced scheme. Insets of the images are defined by the zone in (a) marked out with a white square. 


contain and then reduce the accuracy of the vector cor¬ 
relation calculation m- Another two sets of images are 
collected below the 0 degree and beyond the 180 degree 
positions in order to complete the data sets of the first, 
second, penultimate, and last projections. 

A reference image g(p,x,y,k p ), i.e., without the sam¬ 
ple inserted into the beam, is taken before each 20 of the 
scans. Hence, a total of five reference matrices of size 
2048 x 1850 x 20 could be built by a circular permuta¬ 
tion of the speckle-vector elements and used as reference 
vectors r w (£l n ,x,y)\ 

Y v (QmX,y)= (A), x ,y ,kp) 

{g(( q (C),x,y,k P )} k (2) 

{g(C(E),x,y,k p )} k ) 

where £ is the elementary cyclic permutation over the 
elements of p with order q = mod (n — 1,5). From these 
speckle vectors of length 20, the absorption and dark-field 
signals in addition to the two differential phase gradients 
could be calculated for each Q, n using the XSVT method 
sketched in Fig. [lj (c) and detailed in Ref. [14]. 

In concrete terms, the displacement of the speckle vec¬ 
tors v(Q n ,x,y) induced by the sample refraction is cal¬ 
culated using the Pearson correlation coefficient p\ 

(v x ,v y ) = -axgmaxp(s v (x,y),r v (x + v x ,y + v Y )) (3) 

(v X jV Y ) 


with v(x, y) = v x e x + v y e y and where the variable 
is dropped for clarity. Note that the vectors v are 
calculated with subpixel accuracy [2Tj . 

Next, the differential phase V0 is obtained using the 
formula already mentioned above and valid for the small 
angles: 


„ ,, v 27r w(x. y) 

Vcf>(x,y) = -^f" (4) 

where A is the photon wavelength. We call the wave- 
front gradient VW = v(x, y)/d, which corresponds to the 
angular deflection of the rays caused by refraction. The 
transmission signal T for each projection is calculated by 
taking the ratio of the average speckle- and reference- 
vector intensity elements: 


T(x,y) 


S v (x “h Vxi y “i“ 
< r v (x,y) > 


( 5 ) 


Besides, the dark-field signal Df is calculated by using: 


D f(x,y) 


a(s v (x + v x ,y + v y )) <r v (x,y) > 
<s V (x + v x ,y + v v )> cr(r v (x, y)) 


where a denotes the standard deviation of the vector 
elements. 

Note that this last formula is based only on statistics. 
Hence, in the case of large vectors (r v ,s v ), the approach 
becomes also valid in the case where s v and r v are built 
from images with different illuminations, i.e. with the 
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Absorption (arb. units) 


FIG. 4. The cut in the volume reconstructions of the juniper- 
berry (a) from the phase signal and (b) from the dark-field sig¬ 
nal. (c) An example horizontal transverse slice reconstructed 
from the dark-field signal and the same slice calculated with 
(d) the absolute phase shift signal and (e) the absorption sig¬ 
nal. (f) Cuts along the intensity lines shown in (c,d,e) after 
shift and scaling for comparison purposes. 


membrane located at uncorrelated transverse positions 
between the reference and sample image stacks [27] . 

The phase-imaging method employed here is equivalent 
to x-ray grating interferometry or analyzer-based imag¬ 
ing as the measured signal is the differential phase signal 
V0 [6l |H1 EH]* In these methods, a specific filter is often 
used in the inverse Radon transform to account for the 
specificity of the signal and perform a one-dimensional 
integration in the Fourier space 29, 30]. Here, as the 
XSVT technique provides the two transverse differential 
phase maps, in contrast to the previous techniques, the 
phase images can be recovered by 2D integration by ma¬ 
trix inversion via, for instance, the Cholesky decomposi¬ 
tion. 


A set of raw, directional phase gradient, phase and 
dark-field images is shown in Fig. !(a- h) for the first 
angular projection of the juniper-berry sample. Figure 
[3j (c-f) makes us realize that whilst the traditional XST 
fails largely to track the speckle subsets due the pres¬ 
ence in the sample of high spatial frequency features, 
the XSVT method succeeds perfectly to sense the very 
turbid wave-front. With the proposed interlaced XSVT 
technique, both the vertical and horizontal phase gradi¬ 
ents are correctly recovered for each pixel, which means 
that upon 2D integration, the correct phase images are 
conserved despite the presence of fast-varying structures 
as seen in the inset of Fig. [3j(g). Similarly, the statistics 
offered by the speckle vectors permits a clear sensing of 
the dark-field signal as seen in Fig. [3j (h). 

From the three imaging modalities, 3D reconstruc¬ 
tions are operated using the ESRF PyHST (Python High 
Speed Tomography) [31] implementation with the classic 
filtered back-projection algorithm. Volume reconstruc¬ 
tions from the phase and dark-field signals are shown in 
Fig. [4j(a-b). Transverse reconstruction slices are also 
displayed for comparison in Fig. [4j(c-e). These three sig¬ 
nals do not render the same contrast for numerous sam¬ 
ple features. For instance, the scattering features in the 
outer part of the juniper-berry (noted 1) render a strong 
contrast in the dark-field signal (slice in Fig. [41(c)) while 
a low contrast is observed in the phase slice (Fig. [4j(d)). 
Similarly, within the slice of the phase signal, homoge¬ 
neous features observed within the yellow circle exhibit 
a large phase shift while presenting very low absorption. 
The comparison cuts in Fig. [4| (f) provide a different view 
on these results. 

While the total number of images collected (7200) may 
seem large, it could be readily diminished by reducing the 
number of projections in the angular tomography scan. 
This figure could also be compared with holotomography 
for instance, which also requires a minimum of four im¬ 
ages per projection or with grating interferometry, which 
demands phase-stepping scans of four steps per projec¬ 
tion. For this latter technique, although a special scheme 
could slightly reduce the number of images necessary for 
reconstruction [32], the presence of an absorption grating 
is unfavorable when the irradiation dose received by the 
sample is an issue. The presented scheme can also ben¬ 
efit from a smaller ensemble of images per projection by 
using images of further projections in order to maintain 
the length and then the amount of statistics contained in 
the speckle vectors. 


V. A MIXED XST-XSVT APPROACH 

We present in this section a second scheme for phase 
sensing and tomography reconstruction with the aim of 
reducing the number of acquisitions per projection even 
further. The principle of this method consists of a mixed 
scheme between the traditional XST technique mm 
and the optimized XSVT scheme presented in the pre- 
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Reference Image Stack 



Sample Image Stack 



FIG. 5. Mixed XST-XSVT data-processing scheme. Therein, 
the pixel values considered for the construction of the speckle 
vectors are marked out with black and correspond to the loca¬ 
tion of the central pixel. For the construction of the vectors, 
values of the closest surrounding pixel are used in addition to 
the central pixel value of the five images. 


vious section. The experimental protocol consists of op¬ 
erating the data acquisition according to the procedure 
sketched in Fig. [ 2 ] (top), while recording only one image 
per projection angle, i.e. k p = 1. This means that dur¬ 
ing the scan the light is now modulated with five different 
speckle patterns, with a cyclical repetition of every five 
angular projections. 

As done earlier, the five interlaced images centered on 
each considered projection are used to build a vector with 
sufficient statistics to achieve an accurate tracking. Fur¬ 
thermore, we enlarge the vectors by using not only the 
values of the centered pixels of each image but also the 
ones of the eight surrounding pixel values as displayed in 
f^. m This approach permits to generate, in our case, 
vectors containing 35 components: 


SvO ,y) = ( \ i(to n - 2 ,x + e x ,y + e y )\ 

\ K J (ex,e y ) 

<i(£l n ,X + £ x ,y + £ y ) } ( 7 ) 

l ) (,£ x ,£y ) 

^(^n+2> X + £ x , y + £ y ) | ) 

J (£x,£y)/ 

with here (e x ,e y ) £ {—ft, 0,/i} 2 , h being the pixel size. 
Reference stacks are built in the same manner from the 
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FIG. 6. Images of the juniper-berry sample obtained with the 
mixed XST-XSVT approach, (a) Vertical and (b) horizontal 
wave-front gradients, (c) Phase and (d) dark-held images. 

reference images: 

r v (n n ,x,y) = ( \g(C q (A),x + £ x ,y + £ y )\ 

J (s Xi e y ) 

\g{( q (C),x + £ x ,y + £ y )\ 

\g(( q (E),x + £ x ,y + £ y )\ J 

L J (ex,ey)/ 

(8) 

The following step of the processing method is the one 
of XSVT described by Eq. [3][6] applied to these new vec¬ 
tors. 

The calculated differential phase and phase images ob¬ 
tained with the mixed XST-XSVT technique are shown 
in Fig. |6j (a-d) for the juniper sample. Figure [7] shows 
the equivalent calculated maps for the red-currant sam¬ 
ple. Images (e-f, i-j) are extracted with the interlaced 
XSVT technique previously presented while images on 
the right side, labeled (g-h, k-1), are the maps calculated 
with the present mixed XST-XSVT technique. From 
these two sample projection images, although the two 
methods apparently render similar results, a meticulous 
observation shows that the images obtained with the 
mixed XST-XSVT method reveal less detail. This dif¬ 
ference is better emphasized when comparing the insets 
of the wave-front gradient maps. In the juniper-berry 
sample, some features are smeared out with the mixed 
technique, thus preventing a correct three dimensional 
reconstruction. Likewise, the scattering image obtained 
with the mixed method offers less signal-to-noise ratio be¬ 
cause the statistics required for its calculation have less 
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FIG. 7. Images of the red-currant sample, (a) Raw image of the sample at the first angular projection with an inset showing the 
area marked out with a white square, (b) The corresponding calculated absorption image. Vertical and horizontal wave-front 
gradients calculated with: (c-d) the XST technique, (e-f) the interlaced XSVT method, (g-h) the mixed XST-XSVT approach 
presented in this section, (i) Phase and (j) dark-field images extracted with the interlaced XSVT method, (k) Phase and (1) 
dark-field images extracted with the mixed XST-XSVT approach. 


independent data. 

Conversely, the red-currant sample, which generates 
much less scattering than the juniper-berry sample, could 
be easily reconstructed by using the classic filtered-back- 
projection algorithm from the different signals extracted 
with this mixed approach. Reconstruction samples are 
shown in Fig. [5] This example demonstrates the correct¬ 
ness and greater appropriateness of this approach over 
the traditional XST technique which, in this particular 
example, failed to accurately 3D reconstruct the smaller 
features. 

More generally, the results presented in Figs. S© 
[7] and [5] demonstrate a better fidelity of the signal re¬ 
covery with the interlaced XSVT scheme than for the 


other speckle approaches. In comparison to the tradi¬ 
tional XST approach m, the mixed scheme process¬ 
ing is highly beneficial to the robustness of the phase 
sensing. Where fast and sharp variations of the sample 
electronic density generate wave-front features with sizes 
comparable to the speckle, the direct XST of speckle sub¬ 
sets becomes quickly limited and/or the spatial reduction 
drastically reduced. In this regard, the interlaced mixed 
method offers much higher efficiency and reliability albeit 
at the cost of a slight degradation in spatial resolution 
with respect to the XSVT method. 

This second scheme, integrating both multi-image and 
multipixel values together in a single vector, allows the 
statistics necessary for accurate correlation and tracking 
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FIG. 8. Volume rendering of the tomographic reconstructions of th e red-currant sample (a) from the phase signal and (b) from 
the dark-field signal extracted with the XSVT method of Sec. IV (c) Slice reconstruction from the phase signal and (d) from 
the dark-field signal using the same method, (e) Volume rendering of the reconstruction from the phase signal extracted with 
the mixed XST-XSVT method. Slice reconstructions obtained with this same method: (f) absorption signal, (g) phase signal 
and (h) dark-field signal. 


of the speckle displacements to be improved. Neverthe¬ 
less, such artificially enlarged vectors do not improve the 
statistics as much as one would expect by a straightfor¬ 
ward calculation based on the ratio of integrated pixels 
(here nine). The fact that the values of the neighbor¬ 
ing pixels are partially correlated results in a statistics 
gain lower than the number of pixels considered. This 
can be illustrated by comparing the sensitivity of the 
two techniques when measuring the standard deviation 
of the wave-front gradient maps in areas of 100 x 100 
pixels in the absence of a sample in the x-ray beam. For 
the XSVT method of Sec. |IV| the measured standard 
deviation is of ^0.35 /irad, while it is of ^0.75 /irad for 
the mixed XSVT-XST method. 


As no additional images are required within this pro¬ 
cessing scheme compared to the XST technique or prop¬ 
agation based methods, this approach presents a high 
potential for the imaging of samples for which the dose 
is a sensitive issue. 


VI. CONCLUSION 

In conclusion, we demonstrate that the XSVT tech¬ 
nique is suitable for high-resolution absorption, phase 
and dark-field contrast tomography. The several images 
collected for each projection permit us to overcome the 
limited spatial resolution of the XST technique, whilst 
the interlaced scheme and the use of a nonabsorbing 
scattering membrane allows one to keep the dose at a 
moderate level. Furthermore, the alternative processing 
scheme presented in the second part is a very promising 
way of improving the performance of the XST technique 
in terms of robustness and spatial resolution. Overall, 
the XSVT method is expected to find many applications 
at laboratory sources, synchrotron, and also Compton- 
accelerator-based x-ray sources, for samples where the 
dose and the image resolution are a problematic combi¬ 
nation. 
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